Evidence for anisotropic polar nanoregions in relaxor Pb(Mgi/ 3 Nb 2 /3)03: A neutron 
study of the elastic constants and anomalous TA phonon damping 
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We use neutron inelastic scattering to characterize the acoustic phonons in the relaxor 
Pb(Mg 1/ / 3 Nb 2 /3)03 (PMN) and demonstrate the presence of a highly anisotropic damping mecha- 
nism that is directly related to short-range, polar correlations. For a large range of temperatures 
above T c ~ 210 K, where dynamic, short-range, polar correlations are present, acoustic phonons 
propagating along [110] and polarized along [110] (TA2 phonons) are overdamped and softened across 
most of the Brillouin zone. By contrast, acoustic phonons propagating along [100] and polarized 
along [001] (TAi phonons) are overdamped and softened for only a limited range of wavevectors, q. 
The anisotropy and temperature dependence of the acoustic phonon energy linewidth V are directly 
correlated with the neutron elastic diffuse scattering cross section, indicating that polar nanoregions 
are the cause of the anomalous behavior. The damping and softening vanish for q — > 0, i. e. for 
long-wavelength acoustic phonons near the zone center, which supports the notion that the anoma- 
lous damping is a result of the coupling between the relaxational component of the diffuse scattering 
and the harmonic TA phonons. Therefore, these effects are not due to large changes in the elastic 
constants with temperature because the elastic constants correspond to the long- wavelength limit. 
We compare the elastic constants we measure to those from Brillouin scattering experiments and 
to values reported for pure PbTi03. We show that while the values of C44 are quite similar, those 
for Cn and C12 are significantly less in PMN and result in a softening of (C11-C12) over PbTi03. 
The elastic constants also show an increased elastic anisotropy (2C44/(Cn-Ci2)) in PMN versus 
that in PbTiOs. These results are suggestive of an instability to TA2 acoustic fluctuations in PMN 
and other relaxor ferroelectrics. We discuss our results in the context of the current debate over the 
"waterfall" effect and show that they are inconsistent with acoustic-optic phonon coupling or other 
models that invoke the presence of a second, low-energy, optic mode. 

PACS numbers: 74.72.-h, 75.40. Gb 
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I. INTRODUCTION 

Relaxor ferroelectrics have received a great deal of 
attention from the physics and materials science com- 
munities because of their exceptional dielectric prop- 
erties and potential for use as piezoelectric devices^— 
Pb(Mg 1/3 Nb 2/3 )03 (PMN), Pb(Zn 1/3 Nb 2/3 )03 (PZN), 
and Pb(Zri_ a; Ti 2: )0 3 (PZT) are prototypical relaxors 
that display a broad and frequency-dependent peak in the 
dielectric response^^ Single crystals of solid solutions be- 
tween PMN, PZN, or PZT and ferroelectric PbTi0 3 (PT) 
exhibit piezoelectric coefficients and electro-mechanical 
coupling factors that are much larger than those of con- 
ventional piezoelectric materials. Despite many theoret- 
ical and experimental studies of these and other relaxor 
ferroelectrics, neither the relaxor ground state nor the 
relaxor transition are well understood, and there are few 
materials with comparable dielectric properties^ 

The structural properties of PMN and PZN display 



several common features. At a temperature T^, far above 
the critical temperature Tc below which long range fer- 
roelectric order can be induced in both PMN and PZN 
by cooling in a sufficiently strong external electric field, 
evidence of the formation of local regions of polar or- 
der, also known as polar nanoregions, has been obtained 
from refractive index measurements and neutron scatter- 
ing pair distribution function analysis*^ The existence 
of such locally ordered regions implies that short-range, 
polar correlations are present, and these are manifested 
in the form of strong neutron and x-ray diffuse scattering 
intensity located in the vicinity of Bragg peaks (Ref. \Q- 
TTI ). The geometry of this diffuse scattering in reciprocal 
space is illustrated for two different Brillouin zones in 
PMN at room temperature in Fig. [1] (identical results 
have been obtained for PZN as well). In contrast to the 
sharp, resolution-limited, Bragg peaks that character- 
ize long-range ordered structures, the diffuse scattering 
originating from short-range, polar correlations is broad 
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in reciprocal space and forms rods that extend along 
(110)j 12 ' 13 Several models have been proposed to explain 
this reciprocal space structure, including: pancake-shaped 
regions in real space (Ref . [ill Tl2l. and 1 14) . polar domain 
walls oriented along (110) fRef. I la), correlations between 
chemically-ordered regions (Ref. 16), Huang scattering 
(Ref. [l3 and fill ) , and relatively isotropic displacements 
of the lead cations (Ref. 1 1 9n . Obviously no consensus or 
satisfactory description currently exists of the real space 
structure of the short-range, polar correlations that give 
rise to the large and temperature-dependent elastic dif- 
fuse scattering cross section in these relaxor materials. 

The similarity of the reciprocal space geometries of the 
diffuse scattering in PMN and PZN demonstrates that 
they are structurally similar on short length scales. How- 
ever both compounds exhibit similar long-range ordered, 
average cubic crystal structures as well. At high tempera- 
tures both systems possess cubic Pm3m symmetry, while 
at low temperatures an average cubic unit cell is retained, 
at least insofar as neither compound undergoes a bulk 
structural phase transition as would be revealed by the 
splitting of a Bragg peak. This low-temperature, average 
cubic crystal structure has been well established for PMN 
using both x-ray and neutron scattering techniques* 2 ^— 
However the case of PZN is controversial. Variable en- 
ergy x-ray diffraction studies designed to probe the bulk 
and near-surface regions of single crystal PZN observed 
a cubic phase in the bulk and a rhombohedral phase 
in the "skin", which led to the discovery of an anoma- 
lous skin effect in relaxors j^ 3 - This discovery was later 
challenged by other researchers who, based on neutron 
powder diffraction measurements, concluded that PZN 
exhibits uniform rhombohedral ground state structure 
instead.— On this subject the weight of experimental 
evidence seems to favor an average cubic ground state 
structure for bulk PZN because neutron diffraction mea- 
surements on pure PZN are plagued by enormous ex- 
tinction effects and thus effectively probe only the near- 
surface region, and because similar skin effects were sub- 
sequently reported by x-ray diffraction studies of single 
crystal PZN doped with 4.5% PT and 8% PT (PZN- 
iPT) . 25 ' 26 Moreover, various diffraction studies of single 
crystal PMN doped with 10% and 20% PT (PMN-xPT) 
show that while the near-surface regions (probed by x- 
rays) are rhombohedrally distorted^ the interior or bulk 
of these crystals (probed by neutrons) remains metri- 
cally cubic down to low temperatures, thus confirming 
the skin effect in this systemi 28 ' 29 Subsequent neutron- 
based strain experiments found a significant skin effect 
in large single crystals of PMN as well that has since 
been observed with x-rays ; 26 ' 30 " — and local regions of 
polarization have been directly imaged near the surface 
using piezoresponse force microscopy^ Under the appli- 
cation of strong electric fields the bulk unit cell in PMN 
remains cubic, whereas the structure of the near-surface 
region distorts; however the intensity of the diffuse scat- 
tering arising from short-range, polar correlations can be 
suppressed by an external electric field, but only below 



the critical temperature Tc ~ 210 K, which is indicative 
of a more ordered structure in the materiaL 31 ! 34 

In view of these results we believe that the dielectric 
properties of both PMN and PZN can be described in 
terms of just two temperature scales: a high-temperature 
scale Trf, below which static, short-range, polar correla- 
tions first appear, and a low-temperature scale Tc, be- 
low which long-range, ferroelectric correlations can be in- 
duced by cooling in the presence of a sufficiently strong 
electric field. Each of these temperature scales is re- 
flected in the neutron elastic diffuse scattering cross sec- 
tion in PMN when measured as a function of tempera- 
ture and electric field^ as shown schematically in the 
upper panel of Fig. [5J The high-temperature scale Td 
is defined as the temperature at which the elastic dif- 
fuse scattering cross section becomes non-zero, and the 
low-temperature scale Tc is defined as the temperature 
at which the elastic diffuse scattering can be suppressed 
by electric fields and ferroelectric domains are formed. 
As shown by a recent high-resolution spin-echo study of 
the elastic diffuse scattering in PMN^ the dynamics of 
relaxors also reflects these same two temperature scales. 
The onset of diffuse scattering at Td was found to co- 
incide with the appearance of static (at least on giga- 
hertz timescales) polar nanoregions that coexist with dy- 
namic, short-range, polar correlations. Tc was found to 
coincide with the temperature at which all short-range, 
polar correlations are static (within experimental reso- 
lution)? 3 ^ These two temperature scales are also clearly 
seen in measurements of the bulk polarization in La- 
doped PbZn_ x Ti x 3 (PLZT) where T c marks the tem- 
perature at which the field-cooled and zero-field-cooled 
polarizations diverged These two temperature scales 
can be understood in terms of random dipolar fields 
that are introduced through the disorder that is inher- 
ently present on the B-site of lead-based perovskite re- 
laxorSf 3 ^— As a function of Ti content x, it has also been 
shown that Td and Tc become equal at the morphotropic 
phase boundary (MPB), where the relaxor properties are 
suppressed in favor of those associated with a well-defined 
ferroelectric phase. 41 

While the static structure of relaxors displays a range 
of interesting properties, the dynamics associated with 
the relaxor transition have many unusual features as 
well. The conventional picture of a (displacive) phase 
transition from a paraelectric state to a ferroelectric 
state is characterized by a long-wavelength (q = 0) 
soft transverse optic (TO) phonon for which the fre- 
quency approaches a minimum at Tc, as happens in the 
well-known paraelectric-cubic to ferroelectric-tetragonal 
phase transition in PbTi03 i 42 ' 43 The frequency of the 
soft mode flro is directly related to the dielectric con- 
stant e determined from bulk measurements via the 
Lyddane-Sachs- Teller (LST) relationship, which states 
that 1/e oc {MIto) 2 ■ Therefore, neutron inelastic scat- 
tering measurements of the lattice dynamics provide an 
important tool with which to investigate ferroelectricity. 
In relaxors, the first pioneering work on the lattice dy- 
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namics of PMN was conducted by Naberezhnov et al£± 
Experiments performed by Gehring et al. and Wakimoto 
et al. demonstrated the existence of a soft TO mode in 
PMN and tracked its temperature dependence , 45 i 46 The 
soft mode in PMN reaches a minimum energy around 
400 K and then recovers linearly (i. e. ^l^o °t \T — Td\) 
at lower temperatures. As shown in Fig.[5J 400 K closely 
matches the onset temperature of the elastic (i. e. static) 
diffuse scattering cross section in PMN when measured 
on gigahertz timescales using neutron spin-echo and neu- 
tron backscattering methods* 3 * 5 . Motivated by these re- 
sults, Gehring et al. proposed that the upper tempera- 
ture scale Td, more commonly known as the Burns tem- 
perature, should be normalized from the value of 620 K 
determined from the measurements of the refractive in- 
dex in PMN, to 420 ± 20 K, where the soft TO mode 
reaches a minimum and static ferroelectric correlations 
(i. e. polar nanoregions) first appear4£ 

An important anomalous feature found in these studies 
concerns the TO mode measured near the Brillouin zone 
center, which is unusually broad in energy compared to 
that in PbTiOa and even becomes overdamped, i. e. the 
half-width at half-maximum (HWHM) of the phonon en- 
ergy lineshape T exceeds the phonon energy hflxo over 
a range of temperatures and wavevectors; this indicates 
that some mechanism is present that shortens the life- 
time of long-wavelength TO lattice fluctuations. The 
term "waterfall" was coined to describe this effect, and 
originally the broadening was believed to be the result 
of scattering of TO phonons from polar nanoregions 
However, the same effect is seen in ferroelectric PMN- 
60%PT, which exhibits no diffuse scattering and is not 
a relaxori 5 - - This proves that the anomalous TO broad- 
ening cannot be associated with the relaxor phase but 
may be the result of the disorder introduced through the 
heterovalent nature of the B-site cations. Another ex- 
planation for the broadening was expressed in terms of 
coupling between the TO and TA (transverse acoustic) 
modes as is observed in BaTiOs^ more recently this 
same idea has been used to explain the apparent "wa- 
terfall" effect in PbTe^ However a series of other pa- 
pers have reported the existence of a second, quasi-optic 
mode at low-energies in PMN, which would then com- 
plicate the dynamics near the zone center, thus giving 
the appearance of a broadening or an extra neutron in- 
elastic scattering cross section at low energies . 53 i 54 This 
interpretation is particularly attractive as it could recon- 
cile the apparent discrepancies in the dielectric response 
with the frequency of the low-energy TO mode expected 
from the LST relation described above as well as the 
measured deviations from the expected Curie- Weiss be- 
havior. 55 Still other groups claim that the extra spectral 
weight at low energies originates from quasi-elastic scat- 
tering from short-range, ferroelectric ordering . 38 ! 56 

It is quite clear given the many different aforemen- 
tioned studies that the low-energy neutron scattering 
cross section is neither well understood nor well char- 
acterized in relaxors. In this paper, we provide a de- 



tailed study of the Brillouin zone and temperature de- 
pendence of the transverse acoustic phonons in PMN over 
the entire Brillouin zone. We observe highly anisotropic 
acoustic lattice fluctuations that reflect the existence of 
a coupling between the acoustic harmonic modes and 
the relaxational diffuse scattering. We also provide mea- 
surements of the elastic constants obtained on the tera- 
hertz timescale and compare these with those made using 
other techniques and those of the canonical ferroelectric 
PbTiOa. Our results point to a reduced value of Cn and 
C12 in PMN compared to those in PbTiOa. Further, 
marked changes in the elastic anisotropy (2C44/(Cn- 
C12)) and differences in Cn— C12 are reported here for 
PMN versus those in PbTiOa. We discuss these results 
in terms of the models mentioned earlier in order to de- 
velop a more complete description and understanding of 
the low-energy lattice fluctuations in PMN. 



II. EXPERIMENTAL DETAILS 

All of the neutron scattering data presented in this pa- 
per were obtained using the C5 and N5 thermal-neutron, 
triple-axis spectrometers located at Chalk River Labora- 
tories (Chalk River, Ontario, Canada), the SPINS cold- 
neutron, triple-axis and DCS time-of-flight spectrometers 
located at the National Institute of Standards and Tech- 
nology (NIST) Center for Neutron Research (Gaithers- 
burg, Maryland, USA). Measurements made on C5 con- 
centrated on the low-temperature phonons. Inelastic 
data were measured using a variable vertical-focus, py- 
rolytic graphite (PG) monochromator and a flat PG ana- 
lyzer. The (002) Bragg reflections of the monochromator 
and analyzer crystals were used to analyze the incident 
and scattered neutron energies, respectively. For elas- 
tic measurements, which characterize the nuclear Bragg 
peaks, the PG analyzer was replaced by a nearly per- 
fect (mosaic ~ 0.06° FHWM), single- crystal, germanium 
(220) analyzer. A PG filter was inserted into the scat- 
tered beam to remove higher order neutrons for all ther- 
mal neutron measurements. A sapphire filter was placed 
before the monochromator to remove high-energy fast 
neutrons. The inelastic data have been corrected for con- 
tamination of the incident beam monitor as described 
elsewhere (see Ref. and the appendix of Ref. [H). The 
horizontal beam collimations were either 33'-33'-5-29'- 
72' or 12'-33'-S-29'-72' (S = sample). The final energy 
was fixed to either Ef — 14.6 meV or 14.8 meV. Data 
collected on N5 were taken using a fiat PG monochroma- 
tor and analyzer with a PG filter inserted into the scat- 
tered beam to remove higher order neutrons. Horizontal 
beam collimations were set to 33'-26'-5'-24'-open' with a 
final energy of Ef — 14.6 meV. These configurations pro- 
vided an energy resolution at the elastic line of 0.95 ± 
0.10 meV. A cryofurnace was used to control the sample 
temperature on both C5 and N5, and helium exchange 
gas was inserted into the sample space to ensure good 
thermal contact and equilibrium at all temperatures. A 
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FIG. 1. Elastic diffuse scattering intensity contours measured 
near a) Q = (0,0,1) and b) Q = (1,1,0) at T = 300 K on 
SPINS (E f = 4.5 meV). The sample (S-II) was aligned in the 
(HHL) scattering plane. 



thermometer was also placed on the sample stick near 
the sample position in order to verify and monitor the 
sample temperature. 

Data on SPINS were obtained using a variable vertical- 
focus PG monochromator. A flat PG analyzer was in- 
serted into the scattered beam in series with a beryllium 
filter to remove higher order neutrons, and the final en- 
ergy was fixed to Ef — 4.5 meV. The collimations were 
set to guide-80'-S'-80'-open for inelastic measurements 
and guide- lO'-S'-lO'-open for studies of the nuclear Bragg 
peaks. The energy resolution for the inelastic configura- 
tion was 0.28 ± 0.05 meV. The sample was placed in a 
high-temperature displex so that temperatures from 50 K 
to 600 K could be reached. 

Data collected on DCS were taken using a fixed inci- 
dent energy of 6.0 meV, and the time between pulses at 
the sample was set to 9 ms, which effectively eliminates 
frame overlap. Time-of-flight spectra were collected for 
~ 45 different crystal orientations spaced 0.5° apart at 
each temperature using the 325 detectors with an angu- 



lar coverage from 5 — 140°. More details about the DCS 
instrument can be found elsewhere^ The energy resolu- 
tion at the elastic line was measured to be 0.28 ± 0.03 
meV. The sample was mounted in a high-temperature, 
closed-cycle 4 He-refrigerator such that temperatures be- 
tween 50 K to 600 K could be reached. 

The data presented in this paper were obtained from 
two different single crystals of PMN, which were aligned 
in either the (HK0) or (HHL) scattering planes. The 
room-temperature, cubic lattice constant is a — 4.04A, so 
one reciprocal lattice unit (rlu) equals 27r/a = 1.56 A -1 . 
The crystal labeled S-I is 9 cc in volume, and the second 
crystal labeled S-II is 5cc in volume. We found no evi- 
dence of a structural distortion in either sample. Crystal 
S-I contains a very small amount of PbTiOa that was in- 
tentionally added during sample preparation to facilitate 
the growth of large single crystals; the sample prepara- 
tion method has been described previously^ Crystal S-I 
is the same sample as that used in a previous study of 
the low-energy phonons (see Ref . l6ll) . We emphasize that 
the phonons we measured are identical across all samples 
studied. We compare our measurements on these PMN 
crystals to those on PMN-60%PT, which is not a relaxor 
and which undergoes a well-defined, cubic-to-tetragonal 
phase transition just below 550 KiS . We also compare 
our results to published and unpublished measurements 
made on a smaller PMN single crystal that was studied 
in a series of earlier publications ! 46 ' 47 ' 62 



III. ELASTIC SCATTERING: ANISOTROPIC 
DIFFUSE SCATTERING AND SHORT-RANGE, 
POLAR CORRELATIONS 

In this section, we discuss the temperature depen- 
dence of the diffuse scattering cross section as well as 
its dependence on wavevector when projected onto the 
(HHL) scattering plane. Our results confirm previous ex- 
perimental findings that the diffuse scattering is highly 
anisotropic in reciprocal space. We also discuss how the 
two temperature scales defined in the Introduction are 
manifested in the diffuse scattering data and the conse- 
quences for the polar and ferroelectric correlations. 

The temperature and wavevector dependences of the 
neutron diffuse scattering cross section in PMN have been 
reported in several studies in which measurements were 
made mainly in the (HK0) scattering planei 13 ' 63 These 
studies found that the distribution of diffuse scattering 
intensity in reciprocal space is composed of rods that 
are elongated along (110), which agrees with the results 
of the three-dimensional mapping study conducted with 
high-energy x-rays on PZN doped with PbTiOa. 64 The 
reciprocal space anisotropy revealed by these measure- 
ments has been interpreted in terms of an underlying 
two-dimensional, real space structure in which the long 
axis of the diffuse rod reflects short-range, polar correla- 
tions while the shorter, perpendicular axes reflect longer- 
range correlations. Fig. Q] shows data obtained on the 



5 



cold-neutron, triple-axis spectrometer SPINS that illus- 
trate how the neutron diffuse scattering cross section in 
PMN looks when projected onto the (HHL) scattering 
plane. Panel (a) shows constant-intensity contours mea- 
sured near Q — (0, 0, 1) while panel (b) shows contours 
measured near Q = (1, 1,0). In both cases a butterfly- 
shaped pattern is seen that appears to be slightly elon- 
gated roughly along (111). This is consistent with neu- 
tron measurements performed in the (HKO) scattering 
plane as well as those from the three-dimensional x-ray 
mapping study, which observed rods extending along 
(110). The slight elongation along (111) results from 
rods that protrude out of the scattering plane along 
[110], but which nonetheless contribute to the total dif- 
fuse scattering pattern because of the non-zero vertical 
(i. e. out-of-plane) instrumental Q-resolution. These re- 
sults demonstrate that the diffuse scattering cross sec- 
tion is highly anisotropic in reciprocal space. In particu- 
lar, Fig. [TJ (a) shows that the diffuse scattering measured 
near Q = (0,0,1) is slightly narrower along [001] than 
it is along [110] while the opposite is true for the dif- 
fuse scattering measured near Q = (1, 1, 0) (Fig. [T] (&)). 
As mentioned in the Introduction, this reciprocal space 
anisotropy has been interpreted in terms of a correspond- 
ingly anisotropic real space structure composed of local, 
two-dimensional correlations. 

The temperature dependence of the diffuse scattering 
cross section measured at Q = (0.025, 0.025, 1.05) is sum- 
marized in Fig. [2] (a); this value of Q corresponds to the 
peak in the diffuse scattering intensity when the wavevec- 
tor is scanned along (11,11,1.05). These measurements 
were performed using cold neutrons on SPINS with a final 
neutron energy Ef = 4.5 meV, which provides an excel- 
lent energy resolution (FWHM) of 28 E ~ 0.23 meV. The 
data show a well-defined onset of the diffuse scattering 
intensity at Td ~ 420 ± 20 K that matches the onset mea- 
sured using far better energy resolution (2SE ~ 1 /ueV) 
provided by neutron backscattering techniques. — The is- 
sue of energy resolution is extremely relevant here be- 
cause when the diffuse scattering is measured using a 
technique (e. g. x-rays) that provides a substantially 
coarser energy resolution the resulting intensity will nec- 
essarily include an integration over the low-frequency 
(quasielastic) component of the diffuse scattering cross 
section, and this will lead to an artificially higher (in- 
correct) onset temperature^ Given the invariance of 
the onset temperature of the elastic diffuse scattering 
when changing the instrumental energy resolution by 
more than two orders of magnitude, we conclude that 
Td ~ 420 ± 20 K represents an intrinsic and physically- 
meaningful temperature scale in PMN. The same on- 
set temperature was obtained for several different val- 
ues of Q on different backscattering and neutron spin 
echo spectrometers; hence Td it is not wavevector depen- 
dent. We further note that this temperature matches the 
Curie- Weiss temperature derived from linear fits of the 
inverse dielectric susceptibility measured at a frequency 
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FIG. 2. a) Elastic diffuse scattering intensity measured in 
zero field at Q = (0.025,0.025, 1.05) plotted versus tempera- 
ture. The data were measured on SPINS using a final energy 
Ef = 4.5 meV. The sample (S-II) was aligned in the (HHL) 
scattering plane. The dashed line is a schematic represen- 
tation of the field-cooled data presented in Ref. |3l|. b) The 
thermal expansion coefficient is plotted as a function of tem- 
perature, illustrating no observable structural distortion to a 
sensitivity of a ~ 10 -4 . The data were taken on C5 with 
Ef — 14.6 meV using a perfect germanium analyzer for im- 
proved Q-resolution. The sample (S-II) was aligned in the 
(HKO) scattering plane. 



of 100 kHz ^ it also matches the temperature at which 
the soft TO mode energy in PMN reaches a minimum 
valued 

Cooling PMN below Td does not induce a well-defined 
phase transition to a long-range ordered structure; this 
is evident from Fig. [2] b) , where the thermal expansion 
coefficient of the S-II PMN crystal is plotted from just 
below 200 K to 600 K. Near a conventional ferroelectric 
phase transition, an anomaly is observed in the temper- 
ature dependence of the lattice constant like that seen 
at the ferroelectric transition of PbTiOa^ however no 
such anomaly is seen in PMN. In the relaxor PZN, the 
coefficient of thermal expansion exhibits an anomaly of 
a ~ 5 x 10 -4 at Tc — 400 K in the near-surface re- 
gion^ Efforts to locate a structural distortion in PMN 
were carried out using the C5 thermal-neutron and the 
SPINS cold-neutron, triple-axis spectrometers. A con- 
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figuration that provided excellent wavevector resolution 
was obtained on C5 by closely matching the (i-spacing 
of the sample and that of a nearly perfect Germanium 
(220) analyzer crystal. This configuration and the re- 
sulting resolution function are reviewed elsewhere. 57 The 
data shown in Fig. [2] b) provide no sign of a change in 
the dimensions of the unit cell; instead they imply an 
Invar-like behavior of the lattice constant up to 600 K. 
Therefore, even though Td represents a meaningful tem- 
perature scale in PMN that is associated with the onset 
of static, short-range, polar order and the softening of the 
optic mode, it is not associated with a structural phase 
transition to long-range, polar order. 

In the bulk of single crystal PMN and PZN there is no 
evidence for a well-defined structural transition. Several 
powder diffraction studies of pure PMN have seen sub- 
tle indications of a structural distortion at temperatures 
below 200 K, but these are well-modeled by a disordered 
lattice with local atomic shifts . 20 i 21 i 27 The temperature of 
200 K corresponds closely to that at which a sufficiently 
strong electric held applied along [111] is able to suppress 
the diffuse scattering cross section in PMN (see Fig. [2] 
(a)) while simultaneously inducing a structural distor- 
tion in the near-surface region^ This temperature scale 
is denoted as Tq in Fig. [2] (a). Other studies of single 
crystal PMN and PMN substituted with 10% PbTi0 3 
(PMN-10%PT) using neutron scattering have found no 
evidence of any anomaly in the lattice constant at low 
temperatures, but they do observe a significant change 
above - 350 K (Ref. M and S3). These findings are 
consistent with x-ray diffraction studies of ceramic sam- 
ples of PMN-10%PT^, but they differ from the result 
reported here in Fig. [2] b). Strain-scanning experiments 
performed on single crystals using neutrons have found 
evidence for a significant near-surface effect in PMN and 
even large changes of the coefficient of thermal expan- 
sion as a function of depth in the crystal . 26 ' 30 We there- 
fore conclude that no well-defined structural transition 
occurs at Td ~ 420 K, and that the apparent thermal ex- 
pansion at higher temperatures is sample dependent and 
likely linked to the effects of a skin layer in these ma- 
terials. However, the Invar-like response in the bulk at 
temperatures below ~ 350 K is consistent across a series 
of PMN samples and all known studies. 

Based on these data and those presented in previous 
studies, we conclude that the diffuse scattering cross sec- 
tion is highly anisotropic in reciprocal space. This con- 
clusion is consistent with previous measurements that de- 
scribe the diffuse scattering as being composed of rods 
oriented along (110). 64 Elastic diffuse scattering first ap- 
pears below Td ~ 420 K,— a temperature scale that does 
not reflect a phase transition to a long-range ordered 
structure, but rather to the onset of static, short-range, 
polar correlations. A longer-range ferroelectric phase can 
be induced by cooling in the presence of a sufficiently 
strong electric field, but only below Tc ~ 210 K^i 
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FIG. 3. DCS neutron scattering data measured on PMN sam- 
ple S-II in the (HK0) scattering plane near (200) and along 
[010] at T = 300 K. The constant inelastic scattering inten- 
sity contours reveal well-defined TAi and TO modes near the 
zone boundary at (2, ±0.5,0), but overdamped TO modes 
near the zone center. The vertical dashed lines indicate where 
constant-Q scans were made; these are discussed in the next 
figure. 



IV. INELASTIC SCATTERING: ANISOTROPIC 
TRANSVERSE ACOUSTIC PHONON DAMPING 

So far we have reviewed the wavevector and tempera- 
ture dependence of the neutron elastic diffuse scattering 
cross section in PMN, which reflects the presence highly 
anisotropic, static short-range, polar correlations that de- 
velop below Td ~ 420 K. We now present extensive neu- 
tron inelastic scattering measurements of the transverse 
acoustic (TA) phonons that show how the anisotropy of 
these short-range, polar correlations is manifested in the 
damping of acoustic lattice fluctuations. To this end, we 
discuss the TA phonon scattering cross section measured 
in three different Brillouin zones and the temperature 
dependence observed in each case. 

We chose to make our neutron scattering measure- 
ments in the (HK0) scattering plane because they are 
then sensitive to two different TA phonons, denoted by 



TAi and TA2 . This sensitivity is governed by the neutron 
scattering phonon cross section, which is proportional to 
(Q ' £) 2 ; where Q is the total momentum transfer, or 
scattering vector, and £ is the eigenvector of the phonon. 
Measurements made at reduced wavevectors q transverse 
to (100) or (200), e. g. for Q = (l,q, 0), are sensitive 
to TAi acoustic phonons, which propagate along [010] 
and are polarized along [100]. The limiting slope of the 
phonon dispersion associated with this mode as q — » 
is proportional to the C44 elastic constant. Measure- 
ments made at q transverse to (110) or (220), e. g. for 
Q = (1 + q, 1 — q, 0), are sensitive to TA2 phonons, which 
propagate along [110] and are polarized along [110]. The 
limiting slope of this phonon dispersion is proportional 
to |(Cix— C12). A detailed list of the elastic constant de- 
pendence based on scan direction can be found in Ref. [43| 
and Ref. H 

In this section we will first consider the TAi and TA2 
acoustic phonons measured in the (200) and (220) Bril- 
louin zones as these zones exhibit relatively weak neutron 
diffuse scattering cross sections but strong acoustic and 
optic phonon scattering cross sections. We will then ex- 
amine the TAi and TA2 acoustic phonons measured in 
the (100) and (110) Brillouin zones where strong diffuse 
scattering cross sections are present. This fact has been 
demonstrated quantitatively in PMN at room tempera- 
ture by Vakhrushev et al. , who found that the ratio of the 
diffuse scattering intensities measured in the (110) and 
(220) Brillouin zones was 7 d (110)/7 d (220)>63/4 - 16 
and that between the (300) and (200) Brillouin zones 
was / rf (300)/I d (200)>87/4 - 22. 9 - These ratios should be 
viewed as lower bounds because the experimental deter- 
minations of the weak diffuse scattering cross sections in 
the (200) and (220) zones were limited by experimen- 
tal uncertainties. We mention that, although our data 
cover a wide range of q that spans each Brillouin zone, 
our discussion will not include the highly damped, soft 
optic phonons that are observed at various zone bound- 
aries, because they have already been studied in detail in 
Ref. M- 

A. Q = (2, 0, 0) - TAi acoustic phonons in the 
presence of weak diffuse scattering 

While the (200) neutron diffuse scattering cross sec- 
tion is non-zero, as pointed out in Ref. H3, it is more 
than one order of magnitude smaller than those in the 
(100), (110), and (300) Brillouin zones, which is consis- 
tent with the structure factor calculations discussed in 
Ref. [h] and Ref. [HI We therefore begin our discussion of 
the phonons in PMN in the (200) zone, where the acous- 
tic and optic phonon scattering cross sections are strong 
and the diffuse scattering cross section is weak. 

Fig. [3] shows a contour map of the inelastic scatter- 
ing intensity measured at 300 K along [010] that spans 
the entire (200) Brillouin zone and covers the energy 
range from -20meV to OmeV. The negative energy scale 
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FIG. 4. Temperature dependence of TAi and TO phonons 
measured at two different wavevectors on PMN sample S- 
II in the (HK0) scattering plane using the DCS spectrom- 
eter. The constant-Q scans integrate over the regions (2 ± 
0.025, 0.10 ± 0.02, 0) and (2 ± 0.025, 0.20 ± 0.025, 0). The ver- 
tical dashed lines indicate the average position in energy of 
the TAi phonon. 



merely indicates that these data were obtained using 
a neutron energy gain configuration in which energy is 
transferred from the lattice to the neutron during the 
scattering process; this configuration was chosen because 
it provides a broader detector coverage of momentum- 
energy phase space. At this temperature, well-defined 
TAi and TO phonon branches are seen at large q. The 
TAj mode reaches a maximum energy of ~ 6meV at 
the X-point zone boundary (Q = (2, ±0.5,0)), whereas 
the TO mode is observed at higher energies and extends 
to ~ 20meV. However for wavevectors approaching the 
zone center (q = 0) the TO mode becomes overdamped 
while the TAi remains well-defined. This behavior has 
been termed the "waterfall" effect and was discussed in 
the Introduction^ Previous studies have shown that the 
TO phonon broadens and softens near the zone center as 
a function of temperature, reaching a minimum energy 
at T d ~ 420K4£i£ At temperatures below T d the TO 
phonon hardens and the linewidth narrows. 

This effect is more readily seen in Fig. 31 which displays 



constant-Q scans measured at q = 0.1 rlu and q = 0.2 rlu 
at 600 K, 300 K, and 100 K; these two constant-Q scans 
correspond to the two vertical dashed lines shown in 
Fig. [3] All of the data were fit to two uncoupled harmonic 
oscillators in order to describe the TO mode and the TAi 
mode; the excellent fits are indicated by the solid curves. 
For both small (q — 0.1 rlu) and large (q = 0.2 rlu) 
wavevectors, and at all temperatures, a strong, sharp, 
TAi mode is observed. By contrast, the TO mode is 
only well-defined at large wavevector; at q = 0.1 rlu the 
TO mode is overdamped due to the waterfall effect. On 
cooling from 600 K the intensity of both modes decreases 
in line with expectations based on the Bose factor. A 
detailed analysis of the temperature dependence of the 
structure factors and intensities of the TA modes in PMN 
has already been published^ On the basis of that study 
it was concluded that the TA-TO mode coupling in PMN 
is minimal and does not give rise to the waterfall effect. 
The lack of TA-TO mode coupling is demonstrated by the 
vertical lines in Fig. [H which show that the TAi phonon 
energy does not shift with temperature despite the soft- 
ening and hardening of the optic mode. If strong TA-TO 
mode coupling were present in PMN, as it is in KTa03 
(Ref. l7lh , then a large change in the acoustic mode energy 
position should occur with temperature that mimics the 
change in the TO mode frequency. We observe no such 
change in our experiments. We note that this conclusion 
differs from those of some other studies j 51 i 53 and we will 
discuss these differences at the end of this paper. 

Independent of the issue of mode coupling, Fig. @] 
demonstrates that the TAi mode measured near (200) 
is underdamped (i. e. T<htto) at all wavevectors and 
temperatures studied, and that it does not shift signif- 
icantly in energy with temperature. As demonstrated 
in Ref. [6l], the TA phonon intensities are in reasonable 
agreement with harmonic theory. We note that the TA2 
acoustic phonons measured near (220) also exhibit lit- 
tle to no shift in energy with temperature and are also 
underdamped at all temperatures studied, i. e. between 
100 K and 600 K. In both the (200) and (220) Brillouin 
zones the diffuse scattering cross section is small in com- 
parison to those in other zones such as (100) and (110). 
We now move on to discuss the TA phonon cross section 
in these zones, where the diffuse scattering cross section 
is very strong. 

B. Q = (1, 1,0) - TA2 acoustic phonons in the 
presence of strong diffuse scattering 

The white, dashed line in Fig.[5]represents the set of re- 
ciprocal lattice points where constant-Q scans measured 
in the (110) Brillouin zone are sensitive to TA2 phonons; 
as explained earlier, constant-Q scans directly probe TA 2 
phonons when measured along a direction q transverse to 
Q = (1,1,0). This figure also illustrates the relationship 
between these reciprocal lattice points and the elastic 
diffuse scattering cross section measured at 300 K. The 
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FIG. 5. Elastic diffuse scattering intensity contours measured 
near (110) on SPINS (E f = 4.5 meV) in the (HK0) scatter- 
ing plane. The white dashed line represents the direction in 
reciprocal space where the TA2 and TO phonons were mea- 
sured. The open black circles represent typical reciprocal lat- 
tice points where constant-Q scans were measured. 



diffuse scattering intensity contours were measured by 
Hiraka et at. and are taken from Ref. [63]. The white, 
dashed line also corresponds to a ridge along which the 
diffuse scattering cross section is a maximum compared 
to that along all other parallel lines. Hence constant-Q 
scans measured at points on this dashed line will reflect 
TA2 phonons in the presence of maximum diffuse scat- 
tering intensity. 

Reciprocal space maps spanning the (110) Brillouin 
zone show the TA2 phonon dispersion along [110] at 
600 K and 300 K in Fig.[SJ As mentioned before, because 
of a better detector coverage of momentum-energy phase 
space, our data were taken on the neutron energy gain 
side, which is indicated by the negative value of the en- 
ergy transfers. The constant inelastic scattering intensity 
contours in panel a) provide evidence of a distinct and 
well-defined TA 2 phonon branch at 600 K. But the TO 
phonon branch that was seen in the (200) Brillouin zone 
is not visible here because of a very weak structure factor. 
As in the (200) zone, the TA2 mode reaches a maximum 
energy of ~ 6 meV at the zone boundary and appears to 
be well-defined (underdamped) throughout most of the 
Brillouin zone. The same data are shown after cooling 
to 300 K in panel b). At this temperature the contours 
reveal a TA2 spectrum that is very different from that 
at 600 K because the TA2 lineshape is now overdamped 
throughout most of the Brillouin zone. 

To better understand the temperature dependence of 
the TA2 phonon linewidth T in the (110) Brillouin zone, 
we measured constant-Q scans similar to those shown 
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FIG. 6. Constant inelastic scattering intensity contours mea- 
sured on DCS in the (110) Brillouin zone are shown at (a) 
600 K and (b) 300 K after integrating over (1 ±0.1, 1 ± 0.1, 0). 
Underdamped, dispersive TA2 phonons at 600 K become 
heavily damped at all q at 300 K except for q ~ 0. The TO 
mode is not visible in this zone due to a weak structure factor. 
The wide, black streak at E = meV is where the intensity is 
dominated by the incoherent, elastic scattering cross section 
from the sample and the aluminium sample mount. These 
data were obtained on the S-II PMN crystal aligned in the 
(HK0) scattering plane. 



in Fig. |4j These scans are presented in Fig. [3 which 
summarizes the response of the TA2 acoustic phonon 
at temperatures well above, near, and below X^, where 
the elastic diffuse scattering first appears. The solid 
lines in Fig. [7] are fits to a harmonic oscillator line- 
shape that obeys the principle of detailed balanced The 
dashed lines correspond to a Gaussian function of en- 
ergy centered at E = and describe the elastic scat- 
tering cross section, which is a combination of diffuse 
scattering and incoherent scattering from the sample and 
sample mount. Details of the lineshape and the fitting 
have presented previously in Ref. [6l]. At 600 K under- 
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FIG. 7. Constant-Q scans measured at Q — (1.15, 0.85, 0) and 
(1.4,0.6,0) at temperatures ranging from 300-600 K. The solid 
lines are fits to a damped harmonic oscillator as described in 
the text. The data were taken with the C5 thermal-neutron, 
triple-axis spectrometer located at Chalk River using Ef = 
14.6 meV. The PMN S-I sample was aligned in the (HK0) 
scattering plane. 

damped TA2 phonons are observed at Q = (1.15,0.85,0) 
and (1.4,0.6,0), which are consistent with the contours 
shown in Fig. [6] However after cooling to 400 K the TA2 
mode at Q = (1.15,0.85,0) becomes overdamped and 
resembles a "quasi elastic" lineshape in that it overlaps 
with the elastic line at E — meV. A similar lineshape is 
observed at Q = (1.4, 0.6, 0), which lies very near the M- 
point zone boundary. Finally, at 300 K, which is above 
Tc but below Td where the soft TO mode achieves a 
minimum in energy, the TA2 mode at both wavevectors 
is again described by an overdamped lineshape. 

Figure [5] extends our study of the temperature de- 
pendence of the TA2 linewidth to temperatures be- 
low T c - 200K and wavevectors Q = (1.1,0.9,0) 
and (1.25,0.75,0). Both wavevectors exhibit an under- 
damped TA2 lineshape at 600 K and an overdamped line- 
shape at 300 K; however at 100 K, i. e. below Tc, the TA2 
phonon recovers an underdamped lineshape. This low- 
temperature recovery is more evident at low wavevector: 
at Q — (1.1,0.9,0) the TA 2 phonon lineshape is quite 
sharp and the energy has hardened to a value comparable 
with that measured at 600 K; by contrast the lineshape 
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FIG. 8. Constant-Q scans measured at Q — (1.1,0.9,0) and 
(1.25,0.75,0) at 600 K, 300 K, and 100 K. The solid lines rep- 
resent fits to a damped harmonic oscillator as described in 
the text. The data were taken with the C5 thermal-neutron, 
triple-axis spectrometer (Chalk River) using Ef = 14.6 meV. 
The PMN S-I sample was aligned in the (HK0) scattering 
plane. 



at Q=(l-25, 0.75,0) displays a more subtle recovery inas- 
much as it still resembles an overdamped form. We point 
out that at 100 K (see Fig. [5]) the elastic diffuse scattering 
cross section is maximum. 

These results provide overwhelming evidence of strong 
coupling between the TA 2 mode and the diffuse scat- 
tering centered at E = 0. The fact that strong, 
temperature-dependent shifts in the TA 2 phonon energy 
and changes in the TA 2 linewidth are observed in the 
(110) zone, where the TO phonon cross section is weak, 
but not in the (220) zone, where the TO phonon cross sec- 
tion is strong, proves that these effects are not due to cou- 
pling between acoustic and optic modes. Instead, these 
effects correlate directly with the strength and temper- 
ature dependence of the diffuse scattering cross section. 
This is further corroborated by measurements of the TA 2 
mode in PZN-4.5PT which demonstrated a recovery in 
the acoustic mode when an electric field suppressed the 
diffuse scattering cross section^ It has been shown in 
a series of studies describing the central peak in SrTi03 



(Ref . [73h . and in the appendix of Ref. HT|, that a dynamic 
component of the diffuse scattering can couple to a low- 
energy, harmonic mode to produce an overdamped line- 
shape. Theoretical studies have also discussed this pos- 
sibility^ Recent studies of PMN-30%PT have applied 
this model to try to understand the heavily-damped line- 
shapes observed for the acoustic mode^ Studies of the 
diffuse scattering using neutron spin-echo (NSE) tech- 
niques, which provide a considerable dynamic range in 
Fourier time, have documented the presence of such a 
dynamic component over a wide range in temperatures, 
which becomes static below Tc ~ 210 K. 25, Broad-band 
dielectric measurements suggest that the dynamics may 
even persist to much lower frequencies than are probed 
by the NSE technique. 76 ' 77 To further support the con- 
jecture that the TA phonon damping documented here is 
the result of a coupling between the acoustic and diffuse 
scattering cross sections, we now present measurements 
of the TAi phonons in the (100) Brillouin zone. 



C. Q=(1,0,0) - TAi acoustic phonons in the 
presence of strong diffuse scattering 

The grey, dashed line in Fig. [9] represents the set of re- 
ciprocal lattice points where constant-Q scans measured 
in the (100) Brillouin zone are sensitive to TAi phonons; 
as explained earlier, constant-Q scans directly probe TAi 
phonons when measured along a direction q transverse to 
Q = (1, 0, 0). This figure also illustrates the relationship 
between these reciprocal lattice points and the elastic 
diffuse scattering cross section measured at 300 K. These 
diffuse scattering intensity contours were also measured 
by Hiraka et al. and are again taken from Ref. [63|. How- 
ever, unlike the white, dashed line shown in Fig. [5J the 
grey, dashed line here does not correspond to a ridge 
along which the diffuse scattering cross section is a max- 
imum. Thus the reciprocal lattice points used to study 
the TAi phonon mode in the (100) Brillouin zone are not 
located on ridges of maximum diffuse scattering, as was 
the case for the TA 2 phonons in the (110) Brillouin zone; 
instead they are located between them where the diffuse 
scattering cross section is nonetheless strong. 

We briefly note here that Fig. [3] shows diffuse scatter- 
ing intensity contours measured in the (HK0) scattering 
plane whereas the those in Fig. Q] were measured in the 
(HHL) plane. It can be readily seen that the elastic dif- 
fuse scattering cross-section is composed of rods that are 
extended along the [110] and [TlO] directions. 

Fig. [TU] displays reciprocal space maps, analogous to 
those in Fig. [51 that span the (100) Brillouin zone and 
show the TAi phonon dispersion along [100] at 600 K and 
300 K. The constant inelastic scattering intensity con- 
tours in panel (a) reveal a distinct and well-defined TAi 
phonon branch at 600 K. The same contours in panel (b) 
show that the TAi phonon dispersion does not change 
dramatically and remains well-defined at large q near the 
X-point zone boundary (Q = (0.5,1,0)). This behavior 
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FIG. 9. Elastic diffuse scattering intensity contours measured 
near (100) on SPfNS (E f = 4.5 meV) in the (HKO) scatter- 
ing plane. The grey, dashed line represents the direction in 
reciprocal space where the TAi and TO phonons were mea- 
sured. The open black circles represent typical reciprocal lat- 
tice points where constant-Q scans were measured. 



contrasts with that of the TA2 phonon, which becomes 
overdamped at all q at 300 K except very close to the 
zone center. The data in panel (b) also show that the 
energy of the TAi mode does not shift significantly be- 
tween 600 K and 300 K at the zone boundary. There is, 
however, a strong indication that some extra broadening 
of the TAi phonon is present for wavevectors near the 
zone center at 300 K. 

To explore the wavevector and temperature depen- 
dence of the damping of the TAi phonon in more de- 
tail, we examine the series of constant-Q scans shown in 
Fig. HU At 600 K, which is above both T c and T d , well- 
defined and underdamped TAi phonons are observed. At 
400 K, which is near Td = 420 K, i. e. the temperature 
where elastic diffuse scattering first appears and where 
the TO phonon energy reaches a minimum value, the 
TAi phonons at small wavevector (q = 0.14rlu) become 
overdamped, which is consistent with the behavior sug- 
gested by the contours in Fig. [TD] The TAi phonons at 
larger wavevectors (q = 0.3 rlu) display some evidence 
of a broadened lincwidth, however they remain under- 
damped. This is a key point of this work: the damping 
of the TAi phonon measured at Q = (1,0.3,0) is not 
nearly as dramatic as that observed for the TA2 mode 
at a comparable wavevector measured along [110] in the 
(110) Brillouin zone (see the constant-Q scans measured 
at Q = (1.25,0.75,0) in Fig. IS). At temperatures below 
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FIG. 10. Constant inelastic scattering intensity contours mea- 
sured on DCS in the (100) Brillouin zone are shown at (a) 
600 K and (b) 300 K after integrating over (0,1 ± 0.15,0). 
Underdamped, dispersive TAi phonons at 600 K remain well- 
defined near the zone boundary at 300 K, but become heavily 
damped for intermediate wavevectors. The wide, black streak 
at E = is where the intensity is dominated by the incoher- 
ent, elastic scattering cross section from the sample and the 
aluminium sample mount. These data were obtained on the 
S-II PMN crystal aligned in the (HKO) scattering plane. 

Tc the TAi lineshape recovers, but the phonon is still 
heavily damped at small q. 

The behavior of the TAi phonons measured in the 
(100) Brillouin zone is intermediate to that of TA2 
phonons measured in the (110) Brillouin zone and to that 
of the TA phonons measured in the (200) and (220) Bril- 
louin zones. Both TAi phonons in the (100) Brillouin 
zone and TA 2 phonons in the (110) Brillouin zone be- 
come overdamped on cooling to temperatures near T c i 
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FIG. 11. Constant-Q scans measured at Q — (1.0, 0.14, 0) and 
(1.0,0.3,0) at 600 K, 400 K, 200 K,and 100 K. The solid lines 
represent fits to a damped harmonic oscillator as described 
in the text. The data were taken with the C5 triple-axis 
spectrometer located at Chalk River Laboratory using Ef = 
14.6 meV. The PMN S-I sample was aligned in the (HK0) 
scattering plane. 



(but above Tq) for a range of wavevectors between the 
zone center and zone boundary. However for the TA 2 
phonons this damping extends all the way to the zone 
boundary; for the TAi phonons it does not. The data 
presented in Fig. [TO] show that TAi phonons measured at 
and close to the zone boundary Q — (1, ±0.5,0) remain 
underdamped at 300 K (i. e. below T^). At the other ex- 
treme is the behavior of TA phonons in the (200) and 
(220) Brillouin zones, for which overdamped lineshapes 
are never observed at any temperature or wavevector. 

The low-q damping of the TAi phonons observed in 
the (100) Brillouin zone is consistent with the notion of 
a coupling to the diffuse scattering centered at E = 0, 
and Fig.|9]demonstrates that a significant diffuse scatter- 
ing cross section is present near (100) at small wavevec- 
tors measured along [010]. While we have examined the 
TAi phonons for wavevectors q > 0.14 rlu, it is impor- 
tant to investigate TAi phonons at wavevectors in the 
long-wavelength (q — > 0) limit. It was shown in Ref. 
that TA2 phonons measured in the (110) Brillouin zone 
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FIG. 12. Constant-Q scans measured at Q = (1.0, 1.0, —0.05) 
and (1.0,1.0,-0.15) at 600K, 300K, and 100K. The solid 
lines represent fits to a damped harmonic oscillator plus a 
Lorentzian function centered at E = 0. The data were taken 
with the SPINS cold-neutron, triple-axis spectrometer located 
at NIST using E s = 4.5 meV. The PMN S-II sample was 
aligned in the (HHL) scattering plane. 



at q = 0.10 rlu exhibit a strong damping that persists 
to the Brillouin zone boundary; this finding is consistent 
with the data in Fig.[5J However TA2 phonons measured 
at much smaller q — 0.035 rlu do not show the dramatic 
changes in linewidth and energy that are seen at larger 
q. We observe a similar behavior for the long-wavelength 
TAi phonons measured in the (100) Brillouin zone. Data 
taken with good energy and g-resolution are presented in 
Fig-Hl These data were measured in the (HHL) scatter- 
ing plane near the Q = (1, 1, 0) and display TAi phonons 
at q = 0.05 rlu and 0.15 rlu that can be compared to those 
presented in Fig.[TlJ While a strongly damped TAi line- 
shape is confirmed for q — 0.15 rlu at 300 K, which recov- 
ers upon either heating to 600 K or cooling to 100 K, the 
TAi lineshape at the smaller wavevector of q = 0.05 rlu 
shows comparatively little change in linewidth and al- 
most no shift in energy. We therefore conclude that both 
TAi and TA2 phonons are insensitive to the formation 
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FIG. 13. The temperature dependence of the linewidth F 
(upper panels) and energy Qo (lower panels) of the TAi and 
TA2 phonons was determined by fitting a damped harmonic 
oscillator lineshape to constant-Q scans measured at Q = 
(1.0,0.21,0) and (1.15,0.85,0), respectively. The horizontal 
dashed lines in the two upper panels represent the average TA 
phonon energy and thus indicate when the phonon becomes 
overdamped (r > Qo). The error bars are derived from the 
least-squares fit to the data. 



of polar nanoregions in the long-wavelength limit. The 
essential difference is that the TA2 phonons measured 
in the (110) Brillouin zone are strongly damped all of 
the way to the Brillouin zone boundary whereas the TAi 
phonons are strongly damped only for an intermediate 
range of momentum transfer q. 



V. COMPARISON BETWEEN DIFFERENT 
BRILLOUIN ZONES 

The results from the previous section can be summa- 
rized in the following two figures in which the temper- 
ature dependences of the linewidth (r) and energy po- 
sition (CIq) are plotted for TAi and TA2 phonons mea- 
sured at the same \q\ = V0.15 2 + 0.15 2 ~ 0.21 rlu over 
a large temperature range. Figure [13] presents the line- 
shape analysis of the TA phonons measured in the (100) 
and (110) Brillouin zones, where the diffuse scattering 



cross section is strong, while Fig. [14] presents that for TA 
phonons measured in the (200) and (220) Brillouin zones, 
where the diffuse scattering cross section is weak. 

We start with the (100) and (110) Brillouin zones. The 
upper panels of Fig. [13] show how the linewidth of the 
TA phonons varies from 100 K to 600 K, and the dashed 
lines represent the average TA phonon energy position 
over the same temperature range. These lines there- 
fore represent the threshold above which the TA phonon 
lineshape changes from being underdamped (r<?iOo) to 
overdamped (T > 7if2o)- From this one can immediately 
see that the TA X phonon at this wavevector, although 
heavily damped, remains underdamped at all temper- 
atures. By contrast, the TA2 phonon at the same \q\ 
becomes overdamped between 200 K and 400 K. We em- 
phasize that this does not imply that the TAi phonon is 
never overdamped; on the contrary, the constant-Q scan 
measured at lower q (1, 0.14, 0) shown in Fig. [TT1 reveals a 
clearly overdamped TAi lineshape at 400 K. The point of 
this analysis is to show that for a given temperature the 
overdamping of the TA2 phonon persists to larger q than 
does that for the TAi mode. Hence the damping mech- 
anism affecting the TA modes in PMN is anisotropic in 
reciprocal space, and this anisotropy correlates with that 
of the diffuse scattering cross section. 

The lower two panels of Fig. [13] show that the TA 
phonon energy positions shift significantly with tempera- 
ture. Of particular note is that the TA mode energies are 
both minimal at a temperature that coincides well with 
Trf ~ 420 K (where elastic diffuse scattering first appears, 
see Fig. [5] and the TO phonon energy is lowest). However 
these data also reveal a key difference in that the TA2 
phonon energy position is significantly lower than that 
for the TAi phonon over the entire temperature range. 

Having demonstrated that the TA phonons in PMN 
are strongly coupled to the diffuse scattering cross sec- 
tion, we revisit the (200) and (220) Brillouin zones where 
the elastic diffuse scattering cross section is weak, but 
non-zero. To this end we show in Fig. [TT] the temper- 
ature dependence of the linewidth (r) and energy po- 
sitions (Oo) of the same phonons studied in Fig. 1131 
but measured in the (200) and (220) Brillouin zones. 
Fig. [14] also combines previously published and unpub- 
lished data measured on other PMN single crystals using 
the BT9 thermal-neutron, triple-axis spectrometer, lo- 
cated at NIST, with data taken on the larger S-I and 
S-II PMN crystals studied here. The agreement between 
these data sets prove that there is no sample dependence 
to the results discussed in this paper. They also demon- 
strate that both the TAi and TA2 phonons measured in 
these zones are strongly damped. However, unlike the 
behavior observed in the (100) and (110) Brillouin zones, 
these TA phonons are never overdamped. Moreover, the 
TA phonon energy positions display conventional behav- 
ior in that they slowly increase with decreasing temper- 
ature. We also observe a key difference wherein the TAi 
phonon linewidth narrows markedly below Td while the 
TA2 phonon remains strongly damped to lowest tern- 
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The diffuse scattering in PMN has recently been shown 
to have an inelastic component that was directly observed 
in neutron spin-echo measurements^ and that was pro- 
posed on the basis of cold- neutron scattering experiments 
in Rcf. l56l In these experiments it was found that by fit- 
ting the spectrum the inelastic linewidth of the diffuse 
scattering cross section varied as q 2 . The energy scale 
was found to be comparable to the TA phonon energy, 
thus providing a strong channel for coupling and renor- 
malization of the phonon energy. The diffuse scattering 
cross section along the direction of propagation of TAi 
phonons in the (100) Brillouin zone is not as large as 
that for TA2 phonons in the (110) Brillouin zone. Since 
the energy scale and intensity of the diffuse scattering 
along (100) are presumably much weaker, it seems log- 
ical that it would only provide a coupling over a small 
range of wavevectors where the energy scale matches that 
of the corresponding TAi phonon. The apparent lack of 
coupling at the smallest wavevectors, as illustrated in 
Fig. Q21 would then result from the energy scale of the 
diffuse scattering being too small to strongly couple to 
the TA phonon. This physical model of coupling is con- 
sistent with the general ideas proposed for mode coupling 
in KTaC>3 (Ref. |7l|) where it was suggested that the cou- 
pling constant scales linearly with the reduced wavevec- 
tor q measured relative to the Brillouin zone center. 



FIG. 14. The temperature dependences of the linewidth T 
(upper panels) and energy Qq (lower panels) of the TAi and 
TA2 phonons measured in the (200) and (220) Brillouin zones, 
respectively. The values are obtained from constant-Q scans 
measured at the exactly same reduced wavevector q as those 
shown in Fig. [13] The error bars are derived from the least- 
squares fit to the data; for the energy positions the error bars 
are equal to the size of the data points. 



peratures. It is interesting to note that the largest TA 
phonon linewidths are roughly comparable to those mea- 
sured in the (100) and (110) Brillouin zones; however the 
TA phonons in the (200) and (220) Brillouin zones never 
become overdamped because the energy positions are not 
affected to the same degree as they are in the (100) and 
(110) Brillouin zones. 

The lower-q TA phonons display strong damping be- 
low Td, and it is only below Tq that the TAi phonons 
recover; by contrast the damping of the TA2 phonons per- 
sists to low temperature. This result is consistent with 
the highly anisotropic correlations proposed in several 
models involving domain walls or short-range, polar cor- 
relations existing along [110] as TA phonons propagating 
along this direction will be most disrupted by the disor- 
der. The Brillouin zone dependence of the TA phonon 
energy position that results in a strongly overdamped 
lineshape in the (110), but not in the (220) Brillouin 
zone comes about from the coupling between the acous- 
tic mode and the diffuse scattering cross section. 



It is important to compare our neutron scattering mea- 
surements of the low-energy acoustic phonons in PMN 
with those using other techniques including Brillouin 
scattering and ultrasound. Ref. 78) presents a summary 
of both Brillouin (measured by the authors) and ultra- 
sound measurements (taken from Ref. Fzoh covering a 
broad range of frequencies spanning from the MHz to 
GHz range. Their data demonstrate a softening of the 
elastic constants at temperatures below 300 K, and the 
magnitude of the softening is frequency dependent. The 
linewidth of the Brillouin peaks has also been reported 
in Ref. l80l and a broadening and recovery have been 
reported for the longitudinal acoustic (LA) [001] mode, 
but no broadening was observed for the transverse mode. 
These results are generally inconsistent with the results 
reported here. We observe significant broadening for the 
TAi and TA2 phonons with the former displaying a re- 
covery at low temperatures. We also observe that the 
anomalies in the energy position become less pronounced 
at smaller wavevectors and lower energies. We remark 
that our data reflects the lattice dynamics of PMN on 
the THz energy scale, which is nearly three orders of 
magnitude larger than that reported by ultrasound or 
Brillouin scattering experiments. We therefore conclude 
that the effects reported at much lower frequencies are 
distinct from those presented here. We also note that 
infrared and dielectric experiments have reported signif- 
icant dynamics down to much lower frequencies than are 
accessible with neutron scattering techniques. 
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VI. ELASTIC CONSTANTS OF PMN 

Based on the slopes of the acoustic phonon dispersions 
near the zone center, we can calculate the elastic con- 
stants using the equations of motions outlined in Ref. |43|. 
For a cubic crystal, the elastic constant matrix is deter- 
mined by three values - Cn, C44, and Ci2- C44 is deter- 
mined uniquely from the slope of the TAi mode and Cn 
is fixed by longitudinal modes polarized along [100]. The 
elastic constant C12 is not uniquely determined, however 
the slope of the TA2 mode is determined by the difference 
(C11-C12), and the slope of the longitudinal mode polar- 
ized along [110] is set by the sum (011+012+2044)^ 
Therefore, C12 can be determined from our measure- 
ments once On and C11-C12 are known, and it can be 
checked for consistency between the two branches. 



Material (Technique) 


Cn 


C12 


C44 


PMN(neutron) 


1.39(7) 


0.43(5) 


0.53(3) 


PMN(Brillouin-Tu et al.) 


1.49 




0.68 


PMN(Brillouin-Ahart et al.) 


1.56 


0.76 


0.685 


PT (neutron- Tomeno et al.) 


2.32 


1.06 


0.72 


PT (Brillouin-Li et al.) 


2.35 


0.97 


0.65 



TABLE I. A summary of the elastic constants (in units of 
10 n /V/m 2 ) measured at 300 K for a series of PMN and PT 
compounds. The Brillouin and neutron scattering results on 
PT were taken from Ref. [80l - [83l . The elastic constants were 
calculated assuming a density p=8.19 g/cm . 



We have determined the elastic constants by measur- 
ing the slopes of the transverse and longitudinal acoustic 
phonons near the (200) and (220) Brillouin zone centers. 
We have chosen these zones because the diffuse scatter- 
ing cross section is very weak and the temperature de- 
pendence of the phonon energy positions exhibits little 
effects from the coupling to the relaxational dynamics 
that are prevalent in the (100) and (110) Brillouin zones. 
To extract an accurate value for the limiting slope for the 
phonon dispersion curves we have restricted ourselves to 
wavevectors less than 15% of the Brillouin zone in order 
to assure linearity in the dispersion. 

The values of the elastic constants for PMN at 300 K 
are listed in Table fl] and are compared to values obtained 
from Brillouin scattering measurements. The neutron 
and Brillouin scattering measurements agree well and the 
discrepancy can be resolved through the fact that neu- 
trons probe the dynamics on the THz timescale where 
Brillouin techniques study the dynamics in a significantly 
lower frequency range. The published values for PbTi03 
are also listed for comparison. While C44 for PMN and 
PbTiC>3 are very similar, both Cn and C12 are reduced 
in PMN. The velocity of the TA 2 mode (set by C u -Ci 2 ) 
is also reduced with C u -Ci 2 =0.88 for PMN and 1.32 
for PbTiC>3. It is interesting to note that near the mor- 
photropic phase boundary (MPB), which occurs when 
~ 33%PT is substituted for PMN, the difference C11-C12 



reduces further and nearly approaches zero, but then in- 
creases for Ti-rich compositions beyond the MPB bound- 
ary^ This suggests the presence of an instability of re- 
laxors to TA2 phonons^ This softening of the acoustic 
branch maybe a contributing factor to the increased cou- 
pling observed between the relaxational and harmonic 
modes in the (100) and (110) Brillouin zones. Another 
difference is seen in the values of the elastic anisotropy, 
defined as 2C 4 4/(Cn-Ci 2 ), which for PMN is 1.43 and for 
PbTiC>3 is 0.75. Again, this is indicative of a softening of 
the TA2 branch and possibly indicative of an instability 
to a homogeneous deformation along (110). 

Fig. 03] summarizes the energy positions of the TAi 
and TA 2 modes as a function of temperature in the (200) 
and (220) Brillouin zones. These data show that a small 
increase occurs with decreasing temperature, but there 
is no concurrent softening of the elastic constants. The 
data are suggestive that, on the THz timescale probed 
by neutrons, there are only relatively modest changes in 
the elastic constants. This finding stands in contrast to 
that in Ref.[78|, which presents Brillouin data measured in 
the MHz range that shows a large softening of the elas- 
tic constants (including C44). We emphasize that our 
experiments are limited by the energy resolution of the 
spectrometer to the THz energy range, and we are not 
able to access the MHz region presented in Ref. [7|| As 
discussed in Ref. |86j, the difference in elastic constants at 
these extremes in frequency maybe due to anharmonic ef- 
fects. Dielectric and infrared measurements have shown 
significant dynamics at low energies, which might be con- 
tributing to the apparent discrepancy between the MHz 
and THz ranged 



VII. COMPARISON BETWEEN THE RELAXOR 
PMN AND THE FERROELECTRIC PMN-60%PT 

As demonstrated in Ref. [8^, the diffuse scattering cross 
section can be tuned through PbTiC>3 (PT) doping, and 
a suppression of the diffuse scattering cross section is ob- 
served for large dopings that lie beyond the MPB. Such 
Ti-rich compositions also show no relaxor properties; in- 
stead they exhibit well-defined ferroelectric phases at low 
temperatures similar to that reported in PbTi03. As 
further evidence that the TA phonon damping reported 
here in PMN is related to the polar nanoregions that are 
associated with the diffuse scattering cross section, we 
present several comparative constant-Q scans of the TA 2 
phonons in PMN-60%PT measured near the (110) Bril- 
louin zone in Fig. QjjJ PMN-60%PT displays no measur- 
able diffuse scattering and displays a well-defined, first- 
order phase transition from a cubic to a tetragonal unit 
cell. Concomitant with this transition the zone center 
TO mode softens and then hardens in a manner simi- 
lar to that measured in pure PbTi0 3 4^ The phonons 
in Fig. [12] are also well-defined and underdamped. The 
temperature of 590 K at which these scans were made was 
chosen because it is just above the structural transition of 
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FIG. 15. The TA 2 phonon in PMN-60%PT measured in 
the (110) Brillouin zone. The data were taken using the 
N5 thermal-neutron, triple-axis spectrometer with the sam- 
ple aligned in the (HKO) scattering plane. 
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FIG. 16. A schematic illustration of the phonon dynamics 
for the TAi and TA2 phonons measured in Brillouin zones 
where the diffuse scattering is strong ((100), (110), and (300)) 
and weak (200). The TO mode in the (110) zone was too 
weak to provide any information on the dampening. While 
the TO mode in the (100) zone was also prohibitively weak, 
the diagram above assumes the lineshape should scale as the 
structure factor. Further discussion of the curves is provided 
in the text. 



~ 550 K and serves to demonstrate that the suppression 
of the diffuse scattering cross section removes the strong 
damping and energy renormalization reported earlier for 
PMN in Fig. [HI These results confirm the fact that the 
large damping and softening of the TA2 mode in pure 
PMN is the result of short-range, polar correlations or 
polar nanoregions. 



VIII. DISCUSSION AND CONCLUSIONS 

We have reported the results of a detailed study of 
the acoustic phonons in a number of Brillouin zones in 
PMN. We have observed that the TA phonons measured 
in the (200) and (220) Brillouin zones are well-defined 
and exhibit a conventional temperature dependence ex- 
pected from harmonic theory in the absence of strong 
TA-TO mode coupling. While the TO phonon cross sec- 
tion in these zones is strong and readily observable, the 
diffuse scattering intensity is weak. The situation in the 
(100) and (110) Brillouin zones is completely different. In 



these zones the diffuse intensity is strong and the optic 
phonon structure factor is weak. While the TA phonons 
in these zones (particularly in the (110) Brillouin zone) 
take on a heavily damped lineshape because of a strongly 
renormalized energy scale, the energy linewidth is similar 
to that observed in the (200) and (220) Brillouin zones. 
Because the lattice dynamics in different Brillouin zones 
is invariant, we attribute the damped lineshape in these 
zones to a strong coupling between the diffuse scattering 
intensity and the TA mode. 

For temperatures T C <T < T^, the TA2 phonons stud- 
ied in the (110) Brillouin zone become overdamped 
over nearly the entire zone. However, \ow-q phonons 
representing long-wavelength excitations remain under- 
damped as expected on the basis of a coupling strength 
that scales with The TAi phonons measured in 
the (100) Brillouin zone display significant damping at 
small wavevectors. These phonons recover below T c . At 
larger wavevectors, however, while some damping is seen, 
the TAi phonons remain underdamped at all tempera- 
tures. The difference in behavior between TAi and TA2 
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phonons is consistent with the idea of coupling between 
the diffuse scattering cross section and TA phonons be- 
cause the former is a maximum when measured along 
[110]. 

The behavior of the low-energy TA and TO phonons 
is qualitatively summarized in Fig. 1161 This figure illus- 
trates the low-energy TA and TO modes and denotes a 
broadening of the mode, indicative of damping, through 
the shading. Based on our experimental measurements, 
we observe a direct correlation between the damping of 
the TA mode and the the diffuse scattering cross section; 
we therefore conclude that there is a strong coupling be- 
tween the two. We have also reported strong damping of 
the TA phonon in the (200) and (220) Brillouin zones, 
representing the anisotropy of the polar correlations. 

Indeed the presence of a strong relaxational component 
to the diffuse scattering may help to explain recent hyper- 
Raman data^ which were used to support the case for 
a second, low-energy TO mode. We note that this study 
never observed a well-defined, underdamped mode, but 
rather a strongly-damped feature that could be inter- 
preted in terms of relaxational dynamics. Thus many of 
the anomalous features resulting from extra intensity at 
low energies in the (300) Brillouin zone could be inter- 
preted in terms of the strong diffuse scattering coupling 
with the harmonic modes as observed in our data, which 
were measured in the absence of strong optic modes. 



A. The Origin of the Waterfall Effect is not due to 
TA-TO mode coupling 

We do not agree with the interpretation of the wa- 
terfall effect, or the temperature dependence of the TA 
and TO modes, in terms of coupling between TA and 
TO modes. Such suggestions have been motivated by 
a comparison with BaTiOs^i However, the lattice dy- 
namics of BaTiOs and PMN are completely different. 
The coupling in BaTiOs along [110] results from the fact 
that the TA and soft TO mode branches cross at small 
wavevectors; thus, because these modes also have the 
same symmetry, they become strongly coupled and spec- 
tral weight is transferred between the modes, which then 
produces large energy shifts42 Such a situation does not 
occur in PMN because the lattice dynamics are relatively 
isotropic in comparison to BaTiOs and no mode cross- 
ing occurs. We have shown that no significant changes 
occur in either the TA phonon energy position or the 
spectral weight, both of which should happen in a strong 
coupling picture. While some small amount of coupling 
must exist owing to the fact that the low-energy TA and 
TO phonons have the same symmetry, we do not believe 
that this weak coupling dominates the physics of PMN 
as previously suggested based on the comparison with 
BaTi0 3 . 

We do see anomalies in the TA phonons measured in 
the (100) and (110) Brillouin zones. But we have proven 
that these anomalies cannot be the result of coupling to 



the TO mode because they are absent in the (200) and 
(220) Brillouin zones, precisely where the TO phonon 
structure factor is strong. By contrast, the TO phonon 
structure factor is weak in the (100) and (110) Brillouin 
zones and this phonon is not easily seen in a constant- Q 
scan. This makes the TA-TO mode coupling idea highly 
improbable because the coupling comes from terms pro- 
portional to XFtaFto, where A is the coupling constant 
and F is the phonon structure factor^ A Brillouin zone 
in which the TO mode cross section is weak will there- 
fore display weak coupling effects between TA and TO 
phonon branches. Interpreting the phonon anomalies we 
have observed in (100) and (110) Brillouin zones in terms 
of TA-TO mode coupling is therefore not self consistent. 

Several studies, including Ref. [sT| and Ref. [6^, have 
proposed TA-TO mode coupling in order to explain 
anomalies observed at small wavevectors in the (300) 
Brillouin zone. But this zone has a very strong diffuse 
scattering cross section (see Ref. [Tol ) and a strong TO 
mode, which softens to low energies at small wavevec- 
tors. It is possible that a coupling may exist between the 
diffuse scattering and the TO mode, and this might pro- 
vide an alternate explanation of some of the anomalous 
features in this zone such as the observation of a differ- 
ent g-dependent TO pho non damping compared to other 
Brillouin zones (Ref. 1511 ) and an apparent enhancement 
of the TA mode intensity (Ref. [62|). It might also pro- 
vide an explanation for the presence of a second "quasi- 
optic" mode that was also observed in the (300) Brillouin 
zonei£2 

We believe a much more likely cause of the broadening 
of the soft TO mode (i. e. the waterfall effect) may be 
the presence of strong disorder induced by random, dipo- 
lar fields. This would explain the broadening observed in 
both pure PMN and PMN-60%PT, which is largely ab- 
sent in the parent compound PbTi03. It is interesting 
to note that disorder has been suggested as the origin of 
the broadening of the soft TO mode in BaTiOs4^ 



B. A second low-energy optic mode as the cause 
for the anisotropic dynamics? 

The broadening of the low-energy transverse optic 
mode has been explained in terms of a second mode based 
on hyper-raman fRef. |54f) and some neutron scattering 
results (Refs. andl9lh. We provide an alternative ex- 
planation in this paper for the low-energy spectral weight 
and we do not believe that our results, particularly in the 
Q=(1,1,0) and (1,0,0) zones, are consistent with this pic- 
ture. The broadening that we observe at low energies 
in these Brillouin zones becomes less pronounced and to 
vanish in the limit q — > 0. This is consistent with cou- 
pling to a relaxational component and inconsistent with 
a zone center feature such as a second optic mode. Fur- 
thermore, we note that the broadening coincides with T~ 
400 K where a static component to the diffuse scatter- 
ing is onset at lower temperatures (Ref. 35), therefore 
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linking the anomalous broadening we observe to the re- 
laxational dynamics characterizing the diffuse scattering. 
We note that we have not pursued an understanding of 
the Q=(3,0,0) Brillouin zone (as discussed in Refs. l53l 
andl9lh in this report as the combination of strong optic, 
acoustic, and diffuse scattering are present and will un- 
doubtably complicate the interpretation of the dynamics. 

Hyper-raman experiments (Ref. |54[ ). however, illus- 
trate the presence of significant spectral weight at small 
energy transfers consistent with a second low-energy op- 
tic mode. The temperature dependence was believed to 
point to a softening at 400 K, consistent with the neutron 
results described in Refs. andl9ll. The interpretation 
in terms of two modes is further corroborated by the split- 
ting of higher energy optic modes probed using Raman 
spectroscopy^ - — The low-energy measured peaks, how- 
ever, are significantly broader than the resolution and are 
over damped with the linewidth being larger or at least 
comparable with the energy position. We therefore be- 
lieve that the data is consistent the results presented here 
in the Q=(2,0,0) and (2,2,0) zones where the zone center 
scattering is interpreted in terms of a single over damped 
mode. The hyper raman results maybe also consistent 
with the relaxational component to the diffuse scatter- 
ing observed with spin-echo techniques. The comparison 
made by Raman techniques to higher energy split modes 
is compelling and neutrons have not yet been successful 
in observing a splitting of these higher energy modes. 



C. Anisotropic polar nanoregions 

Several models of the diffuse scattering have challenged 
the idea that polar nanoregions are highly anisotropic in 
real spaced The results presented here challenges these 
ideas and instead point to a damping mechanism that is 
highly anisotropic given the large differences observed be- 
tween the TAi and TA2 phonons in the (HKO) scattering 
plane. The large degree of anisotropy in the damping im- 
plies some sort of real space structure that is also highly 
anisotropic. Phonons propagating along directions that 
exhbiti weak polar correlations are likely to be more heav- 
ily damped than phonons propagating along directions 
where long range correlations exist. Our results further 
suggest that short-range, polar order is responsible for 
the strong damping of the TA phonons and that the size 
of this damping is tied to the correlation length associ- 
ated with the diffuse scattering. We note that the TA 
phonon damping reported here only reflects anisotropic 
correlations. Thus we are not able to favor the "pancake" 
model over the [110] domain model proposed by differ- 
ent groups. The results here simply point to anisotropic 
polar correlations. 



D. Dynamics, two temperature scales, and the 
high-temperature scale Ta 

Considerable debate and discussion have been de- 
voted to the nature and existence of the well-known 
Burns temperature, which for PMN is purported to occur 
around 620 K. Our studies here suggest that this high- 
temperature scale is a purely dynamic phenomena and 
thus depends on the energy resolution of the experimen- 
tal probe used to measure it. This is consistent with 
results from neutron dynamic pair density function anal- 
ysis, which illustrate that the off-centering of the lead in 
PMN is dynamic at temperatures above ~ 350 — 400 K. 95 
Static, polar regions begin to appear at 420 K where we 
observe a minimum in the soft TO frequency and the ap- 
pearance of static diffuse scattering on GHz timescales 
measured with neutron inelastic scattering. This lower 
temperature is also consistent with Raman studies in 
which a splitting of some modes was observed, although 
this was interpreted as evidence for a third temperature 
scaled A third temperature scale has also been postu- 
lated based lattice constant measurements }2! however, 
as demonstrated here, such evidence is highly suspect 
given that the thermal expansion is sample dependent, 
and possibly results from the effect of a significant skin 
effect. 

We believe that all of the results on PMN can be 
broadly interpreted in terms of only two temperature 
scales. A high temperature scale (Tj ~ 420 K), which 
is characterized by the appearance of static polar cor- 
relations and where a minimum in the soft mode en- 
ergy exists, and a lower temperature scale (Tc ~ 210 K), 
where domains form and a long-range ferroelectric phase 
can be induced through the application of a sufficiently 
strong electric field. We believe the Burns temperature 
is dynamical and should not be interpreted in terms of 
a well-defined phase transition. This idea was discussed 
in Ref. [4^ where the Burns temperature was redefined 
to be T d = 420 ± 20 K for PMN, instead of the much 
higher value of 620 K that had been derived from index 
of refraction measurements. 

The two temperature scales discussed here are natu- 
rally reflected in the dynamics of the TA phonons pre- 
sented above. The high temperature scale (Td) is where 
the broadening of both the TAi and TA2 phonons is a 
maximum and where the energy renormalization due to 
the coupling to the relaxational diffuse scattering is most 
apparent. The lower temperature scale (Tc) marks the 
recovery of the TAi phonon in terms of a decrease in 
the damping in energy. This lower temperature scale 
is not strongly apparent in the dynamics of the TA2 
phonon, however, we propose that this is due to the real 
space structure of the low-temperature ferroelectric do- 
mains. We emphasize that our results only demonstrate 
that the ferroelectric correlations in real space are highly 
anisotropic; this is consistent with the pancake model 
(Ref. HI Oj, and 03) 

and all other models that invoke 
anisotropic correlations (for example Ref. KEl and OH). We 
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are not able to distinguish between these various models 
on the basis of our experimental data. 

E. Summary 

We have demonstrated an anisotropic damping and 
energy rcnormalization of the TA phonons in PMN. 
TA phonons in all zones display strong damping be- 
low Td, where static diffuse scattering is onset, in- 
dicative of short-range, polar order. Only below Tc, 
where the diffuse scattering becomes completely static, 
do TAi phonons propagating along [100] recover a nor- 
mal lineshape. However, the TA2 phonons, which travel 
along [110], do not, therefore revealing the underly- 
ing anisotropic structure of the polar correlations. TA 
phonons measured in zones that have a strong diffuse 
cross section exhibit a strong renormalization in energy, 
which is indicative of a coupling to the relaxational com- 
ponent of the diffuse scattering. This coupling might 
provide an explanation for the extra neutron scattering 



intensity measured in the (300) Brillouin zone that was 
used to support the ideas of TA-TO mode coupling and 
the presence of a second, low-energy, TO mode. We have 
also provided measurements of the three elastic constants 
in the THz energy range. These are all in reasonable 
agreement with published data using Brillouin scatter- 
ing. When compared to pure PbTiOa, our data suggest 
an instability exists in PMN that is related to the TA 2 
acoustic branch. 
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